Turbulence in small-world networks 
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The transition to turbulence via spatiotemporal intermittency is investigated in the context of 
coupled maps defined on small- world networks. The local dynamics is given by the Chate-Manneville 
minimal map previously used in studies of spatiotemporal intermittency in ordered lattices. The 
critical boundary separating laminar and turbulent regimes is calculated on the parameter space of 
the system, given by the coupling strength and the rewiring probability of the network. Windows 
of relaminarization are present in some regions of the parameter space. New features arise in small- 
world networks; for instance, the character of the transition to turbulence changes from second order 
to a first order phase transition at some critical value of the rewiring probability. A linear relation 
characterizing the change in the order of the phase transition is found. The global quantity used as 
order parameter for the transition also exhibits nontrivial collective behavior for some values of the 
parameters. These models may describe several processes occurring in nonuniform media where the 
degree of disorder can be continuously varied through a parameter. 

PACS numbers: 05.45.-a, 02.50.-r 



I. INTRODUCTION 

A general scenario for the appearance of turbulence 
in extended systems is spatiotemporal intermittency, 
i.e., a sustained regime characterized by the coexistence 
of coherent-laminar and disordered-chaotic domains in 
space and time m Some of the most extensive in- 
vestigations of this phenomenon have been based on 
model dynamical systems such as coupled map lattices 
1- i ! !>!■ Th e idea is that the ingredients of such 
models — a discrete space, discrete time system of in- 
teracting elements whose states vary continuously ac- 
cording to specific functions — are sufficient to capture 
much of the phenomenology observed in complex spa- 
tiotemporal processes, in particular some relevant fea- 
tures of spatiotemporal intermittency and turbulence. In 
this respect, coupled map systems can be considered as 
mathematically simpler and computationally more effi- 
cient models than partial differential equations of hydro- 
dynamics M . The transition to turbulence via spatiotem- 
poral intermittency has mainly been investigated in net- 
works of coupled maps where the connectivities between 
elements are defined from deterministic rules that provide 
order to their spatial structure. There are examples of 
such studies in Euclidean lattices^, [|, ||, [|, g| , fractal ge- 
ometries H , and treelike arrays ||7 as well as on globally 
coupled systems These investigations have allowed 
the characterization of the transition to turbulence as a 
critical phenomenon in a range of Euclidean and fractal 
dimensions and on several connection topologies. 

Because of their discrete spatial nature, coupled map 
systems seem especially appropriate for investigating 



physical phenomena occurring in heterogeneous or disor- 
dered media. Recently, intensive and interesting research 
has been performed on the theory and applications of 
small- world networks [O, [l2[ [k| Small- world net- 
works are a class of networks with a high degree of local 
clustering and a small characteristic length between any 
two elements. It has been shown that small-world net- 
works describe many natural and artificial networks [ jl5| . 
By varying a parameter, small-world networks can be 
continuously tuned between ordered, deterministic lat- 
tices and completely random networks. In this article, 
we consider coupled maps defined on small-world net- 
works as spatiotemporal dynamical systems. We study 
the nature of the transition to turbulence and the prop- 
erties of spatiotemporal intermittency on these networks. 
We explore the changes induced in those processes as a 
result of the variation in the connection topology of the 
interactions in the system. 

In Sec. II, a general coupled map lattice model for 
the treatment of small-world networks is presented. The 
transition to turbulence in coupled maps on small-world 
networks is investigated in Sec. III. The site map model 
that we employ is based on the one introduced earlier 
by Chate and Manneville for regular Euclidean lattices 
in one and two dimensions and which captures the 
essential features of the transition to turbulence in ex- 
tended systems. Section IV contains the observations of 
nontrivial collective behavior arising in the system. Con- 
clusions are given in Sec. V. 



II. COUPLED MAPS ON SMALL- WORLD 
NETWORKS 
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There are several ways to construct a small- world net- 
work. In this article we employ the small-world network 
construction algorithm originally proposed by Watts and 
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Strogatz |TT|. We start from a ring of N sites, where each 
site is connected to its k nearest neighbors, k being an 
even number. Then each connection is rewired at ran- 
dom with probability p to any other site of the network, 
to avoid self-connections. After the rewiring process, the 
number of elements coupled to each site (which we call 
neighbors of that site) may vary, but the total number of 
links in the network is constant and equal to Nk/2. It 
is assumed that all links are bidirectional. Although this 
algorithm does not guarantee that the resulting graph 
is connected, we have used only connected ones for our 
calculations. 

The state of each site on the network can be assigned a 
continuous variable, which evolves according to a deter- 
ministic rule depending on its own value and the values of 
its connecting elements. We define a coupled map lattice 
dynamical system on a small-world network as 

x t+1 (i)=f(x t (i)) + ^- £ [f(x t (j))-f(x t (i))}; 

(1) 

where Xt(i) gives the state of the site i (i = 1, 2, ... , N) 
at discrete time t; k(i) is the number of neighbors or 
elements connected to site i; S(i) is the set of neighbors 
of site z; e is a parameter measuring the coupling strength 
between connected sites, and f(x) is a nonlinear function 
describing the local dynamics. 

The above coupled map equations can be generalized 
to include other coupling schemes or continuous time lo- 
cal dynamics. Different spatiotemporal processes can be 
studied on small-world structures by providing appropri- 
ate local dynamics and couplings. 



at the value r = 3 in all the calculations. As in ordered 
networks, the transition to the extended turbulent state 
can be characterized through the average value of the 
instantaneous fraction of turbulent sites F t , a quantity 
that acts as the order parameter for the system ||. We 
have calculated (F) as a function of various parameters 
of the system from a time average of the instantaneous 
turbulent fraction F t , as 

<F) = if>. (3) 

Thus (F) = describes a laminar state and (F) E (0, 1] 
corresponds to a turbulent state of the network. Random 
initial conditions were used in all runs with given param- 
eter values of the system Eq.(0). About 10 4 iterations 
were discarded before taking the time average in Eq.(||) 
and T was typically taken at the value 10 4 . Increasing 
the averaging time T or the network size N does not have 
appreciable effects on the results. It should be noticed 
that a minimum number of initially turbulent sites is al- 
ways required to reach a sustained state of turbulence. 

Figure 1 shows (F) as a function of e for different fixed 
values of the probability p in a typical small- world net- 
work with nearest neighbor number k = 10 and size 
N = 10 4 . The error bars shown on (F) correspond to 
plus and minus the standard deviations A{F) (square 
root of the variance) of the time series of F t at each value 
of the coupling parameter e. We have verified, by doing 
25 realizations of the rewiring process described in Sec. 
I, that fluctuations on (F) and on A(F) due to different 
configurations of the networks are not significant. 



III. TRANSITION TO TURBULENCE IN 
SMALL- WORLD NETWORKS 

The phenomenon of spatiotemporal intermittency in 
extended systems consists of a sustained regime where 
coherent and chaotic domains coexist and evolve in space 
and time. A local map possessing the minimal require- 
ments for observing spatiotemporal intermittency in cou- 
pled map systems is |3) 

fto = /§(!- I 1 " 2a; l)> if* e [0,1] () 
J v ' \ x, if x > 1, v ' 

with r > 2. This map is chaotic for f{x) in [0, 1]. How- 
ever, for f(x) > 1 the iteration is locked on a fixed point. 
The local state can thus be seen as a continuum of sta- 
ble "laminar" fixed points (x > 1) adjacent to a chaotic 
repeller or "turbulent" state (x G [0,1]). 

In ordered, deterministic networks, the turbulent state 
can propagate through the lattice in time for a large 
enough coupling, producing sustained regimes of spa- 
tiotemporal intermittency. Here, we investigate the phe- 
nomenon of transition to turbulence in small-world net- 
works using the local map / (Eq. (||)) in the coupled sys- 
tem described by Eq. (nl). The local parameter is fixed 
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FIG. 1: Mean turbulent fraction (F) as a function of the 
coupling e for small- world networks with k — 10 and N = 10 4 . 
The error bars indicate ±1 standard deviations, (a) p = 0. 
(b) p = 0.17. (c) p = p c = 0.55. (d) p = 0.80. 

Figures l(a)-l(d) show that, for a fixed value of p, there 
is a critical value of the coupling e c at which the transi- 
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tion from a laminar state to turbulence occurs. For small 
values of p (Fig. 1(a)), a regime of rclaminarization of 
the system takes place on an interval of the coupling pa- 
rameter, after a window of turbulence. The relaminar- 
ization gap disappears with increasing p, leaving a dip 
in the (F) curve, as seen in Fig. 1(b). For small val- 
ues of p the transition to turbulence as the coupling is 
varied takes place continuously, similarly to a second or- 
der phase transition. As p is increased, the transition 
becomes progressively steeper until it happens discontin- 
uously, as in a first order phase transition. There exists 
a critical value of the probability p c 0.55 at which the 
character of the transition to turbulence changes from a 
second order phase transition (Figs. 1(a) and 1(b)) to a 
first order phase transition (Figs. 1(c) and 1(d)). 

In Figs. 2(a) and 2 (b) we show the mean turbulent 
fraction (F) and its standard deviation A(F) plotted as 
functions of both p and e. Fig. 2(a) shows that the tran- 
sition to sustained turbulence occurs on a critical curve 
on the parameter plane (p, e). The variation in the nature 
of this transition along this critical curve as p increases 
can clearly be appreciated in Fig. 2(a). Typical statis- 
tical deviations are seen for small values of p and e in 
Fig. 2(b); however this figure reveals very large fluctu- 
ations in the instantaneous turbulent fraction occurring 
for larger values of those parameters, and which can also 
be observed as the "bulbs" in Figs. 1(c) and 1(d). As it 
shall be discussed in the next section, this phenomenon 
is associated to the emergence of nontrivial collective be- 
havior in the system. 




FIG. 2: (a) Mean turbulent fraction {F) as a function of p 
and e. (b) Standard deviation A{F) as a function of p and e. 
Network parameters are k — 10, N = 10 4 . 

We have calculated numerically the critical values of 
the coupling e c for the onset of turbulence in small- 
world networks with fixed k — 10 as a function of their 
rewiring probability p. Figure 3 shows the resulting crit- 
ical boundary e c (p) for the transition to turbulence, as 
well as the phase diagram of the system, on the param- 
eter plane (p, e) . The critical coupling value e c increases 
as the disorder in the network, described by the proba- 
bility p, grows. The critical boundary curve e c (p) corre- 
sponds to a continuous, second order phase transition for 
p < p c — 0.55 and to a discontinuous, first order tran- 
sition for p > p c . Figure 3 indicates with a dotted line 
where the first maximum of the mean turbulent fraction 
(F) occurs on the parameter plane. This line of first 



maxima of (F) crosses the critical boundary separating 
the laminar and the turbulent states of the system at the 
value p = pc, and the character of the phase transition 
changes at this point. The rclaminarization gap is also 
indicated on the plane (p, e). 
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FIG. 3: Phase diagram of the system Eq. Q1J) . The critical 
boundary separating laminar and turbulent regimes is shown 
with a continuous line. Thin line: second order phase transi- 
tion; thick line: first order phase transition. The dotted line 
indicates the locus of the first maxima of {F} on the param- 
eter plane. The regions of nontrivial collective behavior are 
bounded by dashed lines. 

For values of the probability p < p c , where a contin- 
uous transition from a laminar regime to turbulence oc- 
curs, the variation of the order parameter (F) near the 
critical curve can be characterized by a critical exponent 
(3 as (F) ~ (e — tc) 13 ■ For fixed p, the exponent (3 can be 
calculated from a log- log plot of (F) vs. (e — e c ). The 
critical exponent (3 varies continuously with the rewiring 
probability p. In Fig. 4 we show the resulting graphs of j3 
as a function of p for several small-world networks char- 
acterized by their neighbor number k. In each case, the 
dependence of the exponent (3 with p is well accounted 
by the linear relation (3 = h(p c — p), where the slope h 
varies with k. As p increases, the exponent (3 becomes 
smaller and the corresponding phase transition from 1am- 
inarity to turbulence gets more abrupt. The change in 
the character of the transition from second order to first 
order should occur at the value p = p c for which the ex- 
ponent (3 vanishes. Figure 4 shows the extrapolation of 
the straight line corresponding to the small-world net- 
work with k = 10 until its intersection with the p axis, 
predicting a critical value p c = 0.55. This, in fact, is 
the critical value of the probability at which the change 
in the character of the transition to turbulence for this 
network was observed in Figs. 1 (c), 2(a), and 3. Wc 
recall that the measures of the characteristic path length 
for the family of small-world networks with k = 10 drops 
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to small values typically associated to a random network 
when the rewiring probability is about 0.55 jyj. Thus 
the critical value p c = 0.55 found for the emergence of 
the first order phase transition should be related to the 
onset of randomness in the network. In fact, it has re- 
cently been reported that the transition from laminarity 
to turbulence in randomly coupled map networks can oc- 
cur discontinuously fllq|. 



Fig. 5 means that the networks need less disorder for 
achieving a first order phase transition to turbulence if 
their neighbor number k is sufficiently large. This reflects 
the fact that networks with large enough k are already 
closer to global coupling, and therefore they require less 
rewiring to behave statistically similar to globally cou- 
pled systems. 
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FIG. 4: Critical exponent for the second order phase 
transtion as a function of the rewiring probability p. k = 4 
(squares); k = 6 (circles); k = 10 (triangles); k — 20 (pen- 
tagons). 

By extrapolating the different lines, one can get pre- 
dictions of the critical probabilities p c for different values 
of k. We have numerically verified these critical values 
p c for several small-world networks possessing different 
neighbor numbers k. Thus the linear relation that arises 
from Fig. 4 seems to characterize the transition to tur- 
bulence via spatiotemporal intermittency in small-world 
networks. Furthermore, Fig. 4 predicts that no critical 
value p c < 1 can be found in the case k = 4, and thus the 
transition to turbulence in a small-world network with 
neighbor number k — 4 should occur continuously (sec- 
ond order) for any p, and should always possess a critical 
exponent j3 > 0. This prediction was also verified numer- 
ically. 

Figure 5 shows the predicted critical probability values 
p c as a function of k. Note that p c decreases rapidly with 
increasing neighbor number k. The characteristic path 
length becomes smaller with increasing k JlT| and, con- 
sequently, the connectivity of the small-world network 
approaches the all-to-all coupling limit of a globally cou- 
pled system, where the transition to turbulence is always 
a first order phase transition [jl0| . In globally coupled sys- 
tems, the absence of spatial relations excludes the pos- 
sibility of supporting small domains of turbulent maps 
which would be necessary for a continuous transition to 
turbulence. The decrease in the value of p c observed in 



FIG. 5: Critical values of the rewiring probability vs. k 



IV. NONTRIVIAL COLLECTIVE BEHAVIOR 

In contrast to the usually expected statistical behav- 
ior, the large amplitudes of the standard deviations ob- 
served on some regions of the parameter plane (p, e) in 
Fig. 2(a) do not diminish with increasing system size 
or with longer averaging time. These large fluctuations 
of a statistical quantity indicate the presence of collec- 
tive motions of the system. For example, Fig. 6 displays 
the return maps of the instantaneous turbulent fraction 
Ft for two different values of the coupling corresponding 
to the two "bulbs" observed in Fig. 1(d). Figure 6(a) 
shows a quasiperiodic orbit in the dynamics of F t in the 
first bulb, while Fig. 6(b) reveals a collective period-two 
motion occurring in the second bulb of Fig. 1(d). The 
large fluctuations measured by A(F) reflect the ampli- 
tude of the collective oscillations of F t that emerge in 
those regions of parameters. This feature corresponds to 
a phenomenon of nontrivial collective behavior already 
observed in the temporal evolution of the mean field of 
several chaotic extended systems j[7) . The quantity Ft is 
a simpler statistical description than the mean field and, 
in this case, it already manifests a nontrivial collective 
behavior of the system over long times. 

The two regions on the parameter plane where states of 
nontrivial collective behavior appear are indicated within 
the turbulent zone in the phase diagram of Fig. 3. Note 
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FIG. 6: Return maps of the instantaneous turbulent fraction 



for a small-world network with k 
0.80. (a) e = 0.67. (b) e = 0.85. 
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that the emergence of nontrivial collective behavior oc- 
curs for large enough coupling and an appreciable de- 
gree of disorder, which itself implies a small character- 
istic length between any two elements in the network. 
Thus, as both the amount of short cuts and the coupling 
strength between different parts of the network increase, 
the information transfer required for the emergence of 
collective behavior is more likely to occur. 

Nontrivial collective oscillations in the turbulent phase 
of coupled map systems have also been observed in high 
dimensional Euclidean lattices |l8| , fractal lattices Q , as 
well as in globally coupled maps [[10| . Our results for 
small-world networks show that ordered connections are 
not essential for the occurrence of nontrivial collective 
behavior. 



V. CONCLUSIONS 

We have investigated a coupled map model for the 
transition to turbulence via spatiotemporal intermittency 
in small- world networks. Although the local dynamics 
is simple, we expect that the essential properties of the 
transition to turbulence in small-world structures is cap- 
tured by this model. Coupled Chate-Manneville maps 
could be regarded as a crude description for the dynamics 
of an excitable medium. The system of coupled maps on 



small-world networks can also be used to study different 
spatiotemporal dynamical processes on these structures 
by providing appropriate local maps or couplings. 

By varying the rewiring probability in small- world net- 
works, the behavior of the transition to turbulence can be 
studied in the regime between ordered lattices and com- 
pletely random networks. The critical boundary separat- 
ing the laminar and the turbulent regimes was calculated 
on the parameter plane (p, e) of the system. We have 
found that the character of the this transition changes 
progressively from a second order phase transition to a 
first order phase transition as the disorder in the net- 
work, measured by the rewiring probability, is increased. 
The critical value of the rewiring probability for the on- 
set of the first order phase transition was predicted from 
the scaling behavior observed in the critical exponent (3 
for small values of the probability. Additionally, we have 
been able to calculate the critical values of the rewiring 
probability as a function of the number of initial nearest 
neighbor in small- world networks. 

Discontinuous transition to turbulence and nontrivial 
collective behavior in the turbulent regime are character- 
istic features of globally coupled maps |fj(i|l . These same 
collective properties emerge in small-world networks as 
their rewiring probability is increased. Because of the 
ubiquity of small- world networks in nature and in human- 
made structures, we may expect to see nontrivial collec- 
tive behaviors arising in many systems if the appropriate 
values of relevant parameters are reached. We may ex- 
pect that other typical phenomena of globally coupled 
systems, such as the formation of clusters of syncronized 
elements 19 , could also be observed in small- world net- 
works. 
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